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Abstract: We study the phase diagram of black D3 geometry with uniformly distributed 
D-instanton charge using the probe D7 brane. In the presence of uniform D-instanton 
charges, quarks can be confined although gluons are not, because baryon vertices are al- 
lowed due to the net repulsive force on the on the probe D-branes. Since there is no scale 
in the geometry itself apart from the horizon size, there is no Hawking-Page transition. As 
a consequence, the D7 brane embedding can encode the effect of the the finite temperature 
as well as finite baryon density even for low temperature. The probe D-brane embedding, 
however, undergoes a chiral phase transition according to the temperature and density 
parameter. We studied such phase transitions and calculated the constituent quark mass, 
chiral condensation and the binding energy of baryons as function of the density. The 
baryon vertex melting is identified as the quark deconfinement. We draw phase diagram 
according to these transitions. 
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1. Introduction 

One of the difficulty in holographic QCD is to discuss the temperature dependence of 
physical observables in the hadron phase. The problem is the presence of the Hawking-Page 
transition as the deconfinement/confinement transition JI]]: background metric describing 
the low temperature phase does not contain any temperature parameter. Describing the 
hadrons in finite temperature requests having a black hole metric, in which gluons are 
always deconfined. Therefore we are lead to the question whether quark can be confined 
while gluons are deconfined. In fact the degree of freedom of two species are different 
in large N c order and the confinement of gluons and that of quarks should be separate 
phenomena at least in large N c limit. However, in most of the geometric background with 



well defined Hawking temperature, a baryon vertex J2j is not allowed as a finite energy 
solution^, ||], because black D-brane has net attractive gravity. 

To have a phase transition, we need a scale other than the temperature. In QCD, scale 
invariance is broken by the chiral and gluon condensation, and it is also known that both 
of them are induced by the presence of the instanton. There has been huge activity to 
utilize the instanton background in gauge theory and (quark) dynamics in the instanton 
background has been practiced by many the topic directly relevant to our paper is the 
instanton gas/liquid model, where mass generation and the chiral symmetry breaking was 
discussed using the anomaly and fermion zero modes. For a review, see the review article 
by Schafer and Shuryak |5[ and references therein. 

It is also well known that the D-instanton is the gravity dual to the Yang-Mill instanton. 
Therefore, it is very natural to try to use a background dual to the instanton gas. In fact, 
a background which is dual to the unformly distributed D-instanton over D3 was suggested 
by Hong Liu and Tseytlin ||], where they also proposed, naturally, that the D3/D-instanton 
geometry is dual of the N=4 Yang Mill theory with constant gluon condensation with zero 
electric/magnetic gluon field, that is < F^ u >= 0, < TrF^F^ >= cq 7^, where F is field 
strength of the gauge fields. Such D3/D-instanton geometry contains non-trivial dilaton 
giving non-zero value of gluon condensation q, which was is identified precisely as the 
D-instanton density. 

While the chiral symmetry breaking was a natural consequence of the instanton back- 
ground, the confinement, was not directly reachable in the gauge theory instanton gas 
model. One of the interesting consequence of the AdS/CFT is that Liu and Tseytlin ac- 
tually have shown that the background has a confinement with a linear quark-antiquark 
potential. Therefore in the gravity dual, both the chiral condensation and confinement 
turns out to be consequences of the presence of the D-instanton charges, which produce 
effectively net repulsive force on the probe D-branes and the probe Fundamental strings. 

The purpose of this paper is to consider the finite temperature and density of the 
holographic dual of the D instanton gas/liquid. For our purpose, we need to extend the 
geometry of Liu and Tseytlin to the finite temperature. The necessary metric was found 
by Ghoroku et al. |9| in the context completely unrelated to instanton gas model. While 
most of the dilaton gravity solution is singular [|], ||] such that temperature is ill-defined, 
their solution allows a Hawking temperature. The geometry is quasi-confining, namely, it 
has following properties |J: (i) gluons are deconfined: In QCD, one can consider the gluon 
propagator directly and state the (de) confinement. However, in hQCD, one can not look 
at the gluon propagator directly since that is colored object whose treatment is beyond 
the gravity limit. Instead, what one looks at is the spectrum of the dilaton field which is 
dual of the glueball operator. IF that spectrum is well peaked or discrete, this is the signal 
of confinement. If the spectral function does not show any peak or feature, it means they 
are all dis-integrated and gluons are in de-confined state. As a theorem, the presence of 
the horizon does not allow any stable glueball spectrum, because in AdS black hole, quasi- 
normal mode spectrum has large imaginary part, (ii) the quark-anti quark potential is 
Coulomb like for short distance, linear for medium distance and screened for large enough 
distances. The turn-over distance depends on temperature and gluon condensation, (iii) 



Baryon vertex solution is allowed in low enough temperature. 

There is no geometric transition as temperature grows because (i) we do not have any 
compactified direction, (ii) We need to work out the thermodynamics of the bulk theory in 
Einstein frame, in which our geometry is the same as the black D3 brane, (iii) the action 
for the dilaton and axion cancel each other. The difference between the black D3 brane 
with and without D-instantons comes from the probe D-brane dynamics, which should be 
calculated in the string frame where dilaton factor is manifest as an over all factor in the 
metric. It turns out that the Chern-Simons term can not cancel the effect of the dilaton 
unlike the claim of reference |j] , which is the origin of the chiral symmetry breaking in our 
paper. On the other hand, in ||, the flat embedding is prefered for some reason and it was 
attributed to the cancellation of the dilaton effect by the CS term. As we have shown in the 
appendix A, this is not the case and brane embedding can not be flat, which in turn induces 
the chiral symmetry breaking. We attribute such effect to the net repulsive force acting 
on probe D7 created by the D-instanton charges. Such breaking of the chiral symmetry 
is consistent with the gauge theory where zero mode of the fermions in the D-instanton 
background requests chiral symmetry breaking. 

For the treatment of chemical potential, we follow the method first suggested in [12, 13] 



and developed in [14, |15|. There are two types of probe brane embedding: one describes 
quark phase, where the flavor brane touches the black hole horizon and a baryon vertex 
is not allowed. Usually this embedding is called black hole embedding. If there are strings 
connecting the horizon and the D7, the latter is deformed to a spiky brane to touch the 
horizon. The other embedding is the one where D7 brane never touches the horizon. On 
such brane, strings can be attached only if a baryon vertex (BV) is allowed and present so 
that the strings connect the BV to D7. We find that there is a phase transition between two 
embeddings. We also find that there is a transition between the two black hole embeddings. 
The binding energy of the baryon and the melting temperature in this background was 
studied in [10, [ll[] by considering one baryon vertex connected to the boundary at infinity 
by N c strings. We treat the baryonic medium using the method developed in 0] where 
each compact baryon vertex is joined with a D7 probe brane through a funnel and such 
configuration is smeared along the D3 direction. 

One final comment here is about the nomenclature: usually instanton gas is for weakly 
interacting far-separated instanton and liquid is for dense and non-trivial interaction. Here 
we are using the words "instanton liquid" since homogenous distribution would not be 
consistent with the well separated gas configuration. We do not focus on the interaction 
strength at all. 

The rest of the paper is planned as follows. In section 2, the background geometry 
is reviewed and various embedding configuration was studied with zero density and chiral 
condensations are calculated. In section 3, we study how the embedding geometry changes 
as we vary 3 parameters: temperature, baryon/quark density and the gluon condensation. 
We find that there can be two phases within quark phase: one chiral symmetry broken 
and the other chiral symmetry restored. In section 4, by calculating the free energy and 
chemical potential as function of density we study the phase transition between the baryon 
and quark phases as well as the phase transition between the two quark phases. We work 



out the complete phase diagram both in canonical and grand canonical ensemble. In section 
6, we give a conclusion and future directions. In the appendix show that the Chern-Simons 
term can not cancel the effect of the dilaton unlike the reference ||] . 

2. Background geometry and D-brane setup 

Here we will briefly review the background geometry and the probe brane setup within 
it. The geometry is the one that is a finite temperature extension of D3/D-instanton 
background with Euclidean signature given by||. The background has a five- form field 
strength and a axion field which couples to D3 and D-instanton, respectively. The ten 
dimensional supergravity action in Einstein frame is given by [p], llTj 

S = I J JOxJg (r - \ (5d>) 2 + \e™ {d X f - iif 5) ) , (2.1) 

where $ and x denote the dilaton and the axion respectively. If we set x = — e +Xo> the 



dilaton term cancel the axion term in (2.1). Then the action becomes that of the metric 



and five-form. The solution in string frame can be written as 



ds 2 1Q = e*/ 2 



^(f^^ + d^ + j^-^ + R^nt 



e =1 + — log 77^2 , X = ~e + Xo, 



/w = V 1 -(?) 4 ' (2 - 2) 

where R 4 = Airg s N c a' . The constant q denotes the number of D-instanton. From the 
AdS/CFT dictionary, it also represents the vacuum expectation value (vev) of gluon con- 
densation. The dilaton factor diverges at the black hole horizon. However, it does not give 
any effect on the thermodynamics of the bulk theory since the latter should be calculated 
in the Einstein frame where its onshell action is the same as that of the black D3 brane 
solution. The geometry has a regular event horizon and Hawking temperature given by 
T = r T /irR 2 . 

Quark-antiquark potential is derived from the expectation value of a Wilson loop using 
U-shaped fundamental string configuration ||. For a given temperature T, U-shaped string 
touches the black hole horizon and splits into two straight strings at certain separation 
L*(T). This critical distance L* increases as temperature decreases as in the usual black 
hole geometry. The geometry corresponds to deconfined phase from the gluon point of view. 
In usual black D3 brane geometry, qq potential is Coulomb like for small separation and 
flat when separation is larger than L*, where it is completely screened. In this background, 
the potential is Coulomb like at short distances but linearly grows until U-shaped string 
touches the horizon. At zero temperature the potential is linear up to indefinitely large 
distance. This is a similarity to the real QCD. But the gluon is deconfined which is different 
from the real QCD. We call such phenomena as 'quasi-confinment'. Such quasi-confining 
property comes from the presence of q, which is the value of gluon condensation. It is the 



D-instanton number and it should be closely related to the chiral condensation as well as 
the gluon condensation. 



Introducing a dimensionless coordinate £ by -^ 
can be rewritten as 
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, the background geometry 



ds 2 = e*/ 2 



where r and £ are related by 
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To describe the embedding of the probe D7 brane, we decompose R part in (2.3) into 



ds 2 = e^ 2 



-^ {f(r) 2 dt 2 + dx 2 ) + ^ {dp 2 + p 2 nj + dy 2 + y 2 # 2 ) 



(2.5) 



D7 brane spans (£, x, p) direction and wraps 5 3 , and is perpendicular to y and direction. 
We can set (ft = using the 50(2) symmetry in x 8 ,x 9 plane. Then induced metric on D7 
brane becomes 



J 2 */2 

as£, 7 = e ' 



■£ (/(r) 2 di 2 + dx 2 ) + ^ ((1 + y'W + p 2 tt§) 
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(2.6) 



where y' denotes to dy{p)/dp and £ 2 = p 2 + y 2 . 

The general action of probe D7 brane can be written as sum of DBI action of D7 
brane and the Chern-Simons term. The embedding dynamics of the D7 brane gives the 
dependence of the chiral condensation on control parameters like temperature, density, 
quark mass etc. as we will describe below. 



3. Chiral condensation at zero density 

Here we study D7 brane embeddings at zero density but we first consider the zero tem- 
perature and zero density case to study q dependence of the chiral condensation and then 
consider the finite temperature to study chiral phase transition temperature as a function 
of gluon condensation. 

3.1 Zero temperature limit 

In this limit, 10 dimensional geometry (|2.2|) becomes near horizon geometry of D3/D- 



instanton system which preserve half of supersymmetry ||, 

' R 2 



ds 2 = e*/ 2 
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- 2 (dt 2 +dx 2 ) + ^(dr 2 + r 2 dSl 2 ) 
X 
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(3.1) 



Induced metric on D7 brane can be written as 

r 2 ,.., .^ R 2 



OSny = e ' 



'D7 



-^ (dt 2 + df 2 ) + ^ ((l + y /2 )^ 2 + P 2 <ml) 



(3.2) 



where r 2 = p 2 -\- y 2 . The presence of D-instanton implies that of axion. The axion can 
couple to the D7 brane world volume through the Chern-Simons term. Therefore D7 brane 
action can be written as 

Sm = Sdbi + Scs 

= - m J da 8 e~ V-det( 5 + 2vra'F) + p 7 f d 8 a^C (8)ll ... i8 , (3.3) 

where a is world volume coordinates of D7 brane, g is induced metric on D7 brane and 
CVg) is the Hodge dual 8-form gauge potential of axion field which is form. If we fix D7 
brane position along the fixed <f> direction, the Chern-Simons term becomes locally total 
derivative as we proved in appendix [A]. Therefore it can not affect the equation of motion 
hence we can ignore the Chern-Simons term. This observation makes the key difference of 
our D7 brane embedding to the one described in || and it will be the origin of the chiral 
symmetry breaking. 

Now, the DBI action for D7 brane can be written as 

S D7 = -vi I dtdp e V ^Jl + y' 2 , (3.4) 

with tj = PL7V3VL2, and the equation of motion is 



± 4£L + «AW =0 . (3 .5) 

For a generic value of q, analytic solution is not available. So we look for numerical one. 
For q = 0, the trivial embedding y = constant is a solution which is consistent with the 



result in. I IS] 



If we turn on q, the solution deforms. The bare quark mass m q and chiral condensation 



c are encoded in the asymptotic form of embedding function y(p) pG, 21 1: 



Q 

y(j>) = m q + ~2 H ■ ( 3 - 6 ) 

P 2 

The embeddings with fixed m q are drawn in Figure p]. As q increases, probe D7 brane 
in the centeral region bends upward more and more so that we have non-vanishing chiral 
condensation that is an increasing function of q. The value q corresponds to the gluon 
condensation < TrF 2 > in boundary theory and hence plays the role of parameter of scale 
symmetry breaking. The m q dependence of chiral condensation is drawn in Figure |2|(a). 
With non-zero value of q, the value of chiral condensation goes to finite value in {m q — > 0). 
Therefore, the chiral symmetry is broken. The q dependences of the chiral condensation 
for given m q are drawn in Figure |2|(b). The condensation is increasing function of q but 




(a) 



(b) 



Figure 1: (a) D7 brane embeddings for m q = with q — 0, 0.1, 0.3, 0.5, 1, 3, 5, 10 from bottom 
to top. (b) D7 brane embeddings for m q — 5 with q = 0, 0.1, 0.3, 0.5, 1, 3, 5, 10 from bottom to 
top . 



decreasing one of m q . At large quark mass, chiral condensation is linearly in q. It is 



consistent with the expectation from field theory [23, 21, |2 
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(3.7) 
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Figure 2: (a) m q dependence of chiral condensation with fixed q. (b) q dependence of chiral 
condensation with fixed m n . 



3.2 Finite temperature without density 

In this section, we will study D7 brane embedding in black hole geometry without chemical 
potential. We use the black hole metric ( |2.2| ) and induce metric on D7 brane with (|2,6|). 
Then DBI action for D7 brane becomes 



S D7 = T 7 jdtdpV(p,yWl + y' 2 , 



(3.8) 



where tj = ^^V^Qs and V = e p 3 Lu+ ■ w_. The dilaton factor seems to diverge at the 
black hole horizon. However, the dilaton factor e comes with uj- which contains the zero 



~ (r — rr) at the horizon. This zero is enough to kill the the logarithmic divergence of 
the dilaton. In other words, the dilaton's log singularity does not change the qualitative 
behavior of the brane dynamics near the horizon. 

The equation of motion for the DBI action can be written following form 



V 



1 + y 



12 



d\ogV_ , 
dp 



d logV 
dy 



0. 



(3.9) 



This equation of motion is highly non-linear. We can get solution only in numerical way 
with proper boundary condition (BC). In the presence of black hole horizon, embedding 
of probe D7 can be classified as 'Minkokwski embedding' and 'black hole embedding' |^4j . 
For the Minkowski embedding, we impose BC: y(0) = yo and y'(0) = 0. For the black hole 
embedding, BC is determined by regularity condition of equation of motion at the horizon: 



y(pmin) = yo, y'(pmin) = tan0, 



(3.10) 



where 6 is the angle between p axis and probe brane position at the horizon. In usual black 
hole geometry, gravitational attraction of black hole bends the probe brane downward. 
However, non-zero value of q gives net 'repulsive' force on probe D7. Therefore, D7 bends 
upward, q dependences of Minkowski embedding and black hole embedding are drawn in 
Figure ||[ For fixed quark mass and temperature, the bigger is q, the more pushed up is 
the D7 brane. Such effect is common both in Minkowski and black hole embedding. 

£r=0.7,m,=l 



£j.=l,m,=0.1 
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(a) 



(b) 



Figure 3: (a) q dependence of Minkowski embedding for m q = 1 with q — 0, 1, 3, 5, 10 from 
below, (b) q dependence of black hole embedding for m q = 0.1 with q = 0, 1, 3, 5, 10 from below. 



We expect a phase transition between Minkowski and black hole embeddings as we 
increase temperature. Namely, Minkowski embedding in low temperature will change into 
black hole embedding in high temperature. In the presence of q, D7 feels repulsive force, 
therefore we expect that as q increase, phase transition temperature goes up. One may 
expect that the phase transition might be smoother compared with the case of no instanton 
charge. One can extract q dependence of phase transition temperature Q from the free 
energy. The case for m q = 1 is drawn in Figure [|. 

Now, we will discuss about chiral condensation for this system: it was defined as 
the slope of D7 brane at asymptotic region. See (|3.6[) . We first fix the temperature and 
calculate it with different m„'s. For q = 0, D7 brane bends down and its slope at the 
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Figure 4: q dependence of phase transition temperature. 



asymptotic region is positive and therefore the condensation value is negative. See Figure 
|5|(a). If we turn on q and increase it slowly, the brane bends up relatively very small value 
of q since the pushing-up effect of q is very effective. Therefore the sign of the condensation 
flips at small value of q. Also the value of m q where phase transition occur decreases as q 
increases. 
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Figure 5: m q dependence of chiral condensation c for £r = 1- Blue dashed line denotes to chiral 
condensation for Minkowski embedding and red line denotes to black hole embedding. Dotted line 
indicates phase transition between two embeddings. 



As we decrease m q with fixed q, there is phase transition from Minkowski embedding to 
black hole embedding at certain value of m q . See Figure 0. In Figure || (e,f) there is range 
where m q becomes negative. If we calculate free energy, there is phase transition before m q 
becomes negative. Moreover, in Figure 0(f), free energy of Minkowski embedding is always 



smaller than free energy of black hole embedding. It means that in m q — > limit, the value 
of chiral condensation does not vanish, see Figure 0. This is a spontaneous breaking of a 



U(l) symmetry which is analogue of the chiral symmetry discussed |2C, 21]. Although the 
relevant symmetry is rotation in x 8 ,x 9 plane which is not a true chiral symmetry, one can 



develop the Gellman-Oakes-Renner (GOR) relation [20, 21], which is the purpose of having 
the chiral symmetry. Therefore from now on we call this as chiral symmetry breaking and 
we call c as the chiral condensation. 



fr=i,q=i5 



^>^ ~ 




" C 




r 


B 




2 






A 







-0.3 -0.2 -0.1 



(a) 



(b) 



(c) 



Figure 6: (a) Figure 5(f). A, B and C denote three solution which give zero quark mass, (b) 
Probe brane embeddings for each case, (c) Free energy as a function of quark mass. C (Minkowski) 
embedding has minimum energy. 



For a given temperature, the chiral symmetry is broken, if q is large enough. The q 
dependence of chiral symmetry restoring temperature is drawn in Figure |7| When q — > 0, 
chiral symmetry restoring temperature goes to zero. 




Figure 7: q dependence of chiral symmetry restoring temperature. xSB denotes chiral symmetry 
breaking phase and %S to chiral symmetry restored phase 



4. Chiral symmetry at Finite temperature and finite baryon density 

In this section, we discuss D7 brane embedding with finite density using induced metric 
Q2.6| ). Adding density corresponds to turning on U(l) gauge field A t (p) on D7 brane. The 
DBI action of D7 brane can be written 



S D7 = -t 7 I dtdpp^e^^oj^Je^/ 2 — (l + y 2 )-F 2 := i dtdp£ D7 , (4.1) 



10 



where 

r 7 = n 7 V 4 Q 3 , F = 2ita'F tp (4.2) 

and dot denotes the derivative with respect to p. We used A p = gauge. For a fixed charge 
dynamics, we need Legendre transformation of Lagrangian, which we call 'Hamiltonian': 



2 



n D7 = F^- - C D7 = T 7 Je* — (l +■ f)J& + ^ e * w 3 ( 4. 3) 

<9F y ^+ v 

where Q = Q/(27ra'T 7 ), Q is the number of source charges. Notice that near the horizon, 
the last factor is dominated by the dilatonic term and the whole action is reduced to the 
case of Q = 0. Therefore the argument for the regularity near the horizon goes exactly 
the way of the previous section. We can get numerical solution for the equation of motion 
provided we have proper boundary conditions. The stringy objects corresponding to the 
sources on D7 brane is the end points of fundamental strings. 

Unlike usual black hole background, our background permits the presence of baryon 
vertex. Therefore there are two way of attaching fundamental strings on D7 brane. One 
is connecting D7 to the black hole horizon and the other is connecting it to the baryon 
vertices. Two configurations give different boundary conditions for the D7 brane dynamics. 
We need to examine which configuration has lower free energy. 

4.1 Quark phase 

One way to put point electric sources on D7 brane is to add fundamental strings such that 
one end of strings are on D7 and the other end on black hole horizon. It is equivalent to 
add freely moving quark in boundary theory because the fundamental strings can move 
freely on both D7 and horizon. We call it 'quark' phase. Since the tension of D7 brane is 
always smaller than that of fundamental string p5[ , D7 brane pull down to the horizon. 
Regularity at the black hole horizon requests 

y(Pmin) = tan 0, (4.4) 

where 9 is polar angle of the position where D7 brane touch the horizon. As discussed in 
the previous section, the presence of q gives repulsion on probe D7 brane, and it affects 
its embedding. The q dependence of D7 brane embeddings are drawn in Figure || for two 
different density Q. 

From the figure, we can see that as q increases, the repulsion effect on D7 also increases 
in small density Q. However, if Q is not small, the brane embedding is less sensitive to q 
as shown in Figure |8|(b). This is because the charge Q introduces flux whose electric field 
energy increases the tension of the D7. That is, the stiffness due to the flux is dominating 
repulsion due to the q. 

In the absence of q, embedding solution which corresponds to m q = is uniquely 
determined to be the flat embedding, y(p) = for which value of chiral condensation is 
automatically zero. However, in the presence of q, there are two different embeddings for 
the given m q . For m q = 0, one is trivial with zero chiral condensation c and the other has 
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Figure 8: D7 brane embeddings for m q = 1 and £t = 1 with g = 0, 2, 4, 6, 8, 10 from below for 
(a) Q = 0.1, (b) Q = 10 
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Figure 9: (a) Two embedding solutions with same temperature, density and q. (b) Phase boundary 
between chiral symmetry breaking and restored phases with q = 20. 



non-zero chiral condensation as we can see in Figure ||(a). The solution with non-zero c 
exists in low temperature and small density region. In large density region, only y = is the 
solution. Actually we found that there is a phase transition between the two embeddings 
in certain temperature and density. The phase diagram is presented in Figure ||(b). As q 
increase the phase boundary expands toward larger temperature and density region. See 
Figure [H]. This is natural since the chiral symmetry breaking is caused by the effect of 
gluon condensation q, as we have seen in Figure g 1 . 

In Figure 0, D7 embedding for non-zero m q is drawn. In this case, behavior of D7 
brane embedding is similar to black hole embedding with q = 0. Due to the non-zero value 
of y(oo) ~ m q , chiral symmetry is explicitly broken. However, it is known that there is first 
order phase transition in black hole background between two quark phases in small density 



region[14, 15]. This phase transition line finishes at certain density and temperature. At 
this point, the order of phase transition is second. As q increase, the phase phase boundary 
line moves upward in (T,Q) plane. See Figure |ll]. 



1 This phenomenon is very similar when we turn on the magnetic field on probe brane J2q] in black D3 
brane background. But in that case, the magnetic field cannot affect the background geometry and the 
interpretation of magnetic field also different. 
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Figure 10: phase diagram of chiral symmetry restoration within quark phase for q=15, 25 and 50. 
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Figure 11: q dependence of phase transition between two quark phase for fixed m q = 1. These 
phase transition lines end with second order phase transition point. 



There is an issue about thermodynamical instability around first order phase transition 
point in black hole background ||14| , |l5[ : chemical potential decreases as density increases, 
that is, dfi/dQ < 0. But in the presence of q, this instability can be cured as we will 
discuss later. 

4.2 Baryon phase 

In this section, we study a baryon vertex in the D3/D-instanton background. Baryon vertex 
is originally proposed in Q. Near horizon geometry of black D3 brane is AdS§ xS 5 . If 
we wrap spherical D5 brane on S 5 , due to the Chern-Simons interaction between R-R five 
form field strength and D5 brane world volume, U{\) gauge field is induced on D5 brane 
world volume and to cancel these fluxes we need to put N c fundamental strings on D5 
brane. In asymptotic region, this object looks like bound state of -/V c fundamental string 
or quark. We call this object to 'baryon vertex'. 

As discussed in j3j, however, in Schwartzschild type black hole background does not 
allow compact D5 brane as a solultion of equation of motion of DBI action. Therefore, 
the quark phase is only physical in finite temperature system. But, in the D3/D-instanton 
background, compact D5-brane with N c fundamental strings can be formed even in the 
black hole background. flQ|| . By connecting baryon vertex and probe D7 brane y], we can 
discuss thermodynamics of finite density and temperature system. 
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To study properties of baryon vertex, we rewrite the 10 dimensional metric (| 

„2 /^2 



ds 2 = e^ 2 



-^ (f(r) 2 dt 2 + dx 2 ) + B 2 (^ r+ d9 2 + sin 2 9dtl 



(4.5) 



We take (t,9 a ) as world volume coordinates of a compact D5 brane, and turn on the U(l) 
gauge field on it to have F t e 7^ 0. As the ansatz for the embedding of compact D5, we 
assume the 50(5) symmetry so that position of D5 brane and gauge field depend only on 
9 i.e. £ = £(#), A t = A t {9), where 9 measure the polar angle of S 5 from the north pole. 
The induced metric on D5 brane is 

„2 /t/2 

e 

where £' = d£/d9. The DBI action for single D5 brane with N c fundamental string can be 
written as 

S D5 = -n 5 J e- V-det(s + 27ra'F) + /x 5 I A {1) A G {5) 



1 2 $/2 

ds D5 = e 1 



r 
R 2 



^fdt 2 + R 2 ( ^ + 1 ) d9 2 + R 2 sin 2 9dn\ 



(4.6) 



r 5 / dtd9 sin 4 0e* 



bj_ 



where 



' e *Z^(^2 + ^2)_ F 2 +4At 

= / dtd9£f)5, 
t 5 = fi 5 n 4 R 4 r T , F = 2Tra'F te , A t = 2na'A t . 

fr=0.5, f =2 



(4.7) 
(4.8) 




(a) 



(b) 



Figure 12: (a) £0 dependence of D5 brane embeddings with q = 10, £t = 0.5 . Gray disk denotes 
to black hole.(b) q dependence of tip of D5 brane(£ c ). 



After solving the equation of motion for gauge field and substituting it to the La- 
grangian density, we can get 'Hamiltonian' density of D5 brane, 



H D5 = rj^-— (£' 2 + e)jD(9) 2 + sin 8 9, 



(4.9) 
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where 



D{9) 



+ - sin 9 cos ( 



+ sin 9 cos ( 



(4.10) 



Here, we consider all fundamental strings are attached at north pole. For more detail, see 
Appendix ||. 

The equation of motion for ( |4.9| ) depends on two parameters q, £t and two initial 
condiation £(6> = 0) = £0, £,'{0 = 0) = 0. Here we set 9 = is south pole of D5 brane. For 
a given value of expectation value of gluon condensation(g) and temperature(^), we can 
get numerical solution in terms of £o- The set of numerical solutions in £ plane is drawn 
in Figure |i~2|(a). In this figure, a cusp appears at 9 = it on where N c fundamental strings 
are attached. But these closed D5 brane does not exist in whole value of q and £j. When 
the value of q decrease, the position of tip (we denote it to £ c ) increase and q — >• limit, 



the tip of D5 brane goes to infinity, 12(b). It is consistent with the fact that the usual 
Schwartzschild black hole background does not allow baryon vertex. We also check that 
the baryon vertex solution does not exist at high temperature. 

As we discussed in H , we can add D7 brane at the tip of D5 brane with force balance 
condition |2T 



y'(p 






(4.11) 



where £ c and £' c denote to the position and slop of D5 brane at 9 = it. For details, see 
Appendix |y. The full configurations of D7/D5 brane embeddings with certain value of 
q and density are drawn in Figure [b|. We call this phase as 'baryon phase' because in 
this phase the physical object is baryon vertex. From this figure, we can see that in the 
presence of baryon vertex, the slope of probe brane in asymptotic region is always non- 
zero. It means that in baryon phase, chiral symmetry is always broken. We found that the 
embedding corresponding to the baryon phase does not exist at high temperature. 



fr=0.8, q=15, m,=0, 0=5 



&=0.8, q=15,m,=l,Q=5 




/ \ 
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(a) (b) 

Figure 13: D7 brane embeddings for £r = 0.8, Q = 5 with (1) m q = (b) m q = 1 

5. Phase Transition 

As we discussed in previous section, two kinds of embeddings can exist in finite temperature. 
They correspond to two phases: black hole embedding corresponds to the quark phase and 
Minkowski embedding corresponds to the baryon phase. Quark phase can exist whole 
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temperature region while baryon phase can exist only at low temperature region. In our 
model, the baryon phase exists if temperature is low enough regardless of how high is the 
density. However, in low temperature, quark phase is also possible. To determine which 
is the physical phase, we need to compare the free energies of two systems. There are two 
ensembles we can choose, one is canonical ensemble where density is control parameter. 
In canonical ensemble we can determine physical phase by comparing free energy. Since 
we need to determine the configuration in a fixed value of the charge we need Legendre 
transformed action which we call 'Hamiltonian', although it is not the time translation 
generator. The other one is grand canonical ensemble where chemical potential is control 
parameter. In this case we have to calculate grand potential which is value of DBI action. 

5.1 Canonical ensemble 

To determine physical solution in canonical ensemble, we have to compare free energies of 
two phases. Free energy of quark phase is given by integrating Hamiltonian ( f4.3| ) density 
for the embeddings, 



^quark(Q) = T 7 / d P U D7 (Q) , (5.1) 

Jp min quark phase 



where T-Ldi = ~Hd7/t7 (we introduce it for convenience). Here, we regularized free energy 
by subtracting energy of flat D7 brane. 

On the other hand, baryon phase needs compact D5 brane, called baryon vertex oper- 
ator. Since each D5 should have N c quarks, their number is related to the quark number 
by Nb = Q/N c . The total free energy in baryon phase can be obtained by adding energy 
of D5 brane to that of D7 brane: 



-^baryon(Q) = T 7 / dpH D1 {Q) 

dpHmiQ) 



baryon phase Nq 



+ ^-r 5 / <im m 



baryon phase o7T 



+ jLq I dOUu:, 



(5.2) 



By comparing the value of (5.1) and ( |5.2|) , we can determine which phase is physically 
favored. As we discussed in previous section, there are two phases in quark phase, which 
is true also in massless quark case. One is quark phase with broken chiral symmetry and 
the other one is phase with chiral symmetry restored. The density dependences of free 



energy in massless quark case are drawn in Figure 14. In the figure, we plotted F which 
is defined as F = F — a(T)Q to visualize the difference of the free energies of different 
phases. One should notice that all the free energies monotonically increase as function of 
density, although the figure does not show that due to the subtraction of a(T)Q. 

At low temperature, the free energy of baryon phase is always lower than that of quark 



phase for all density region, see Figure 14(a). It means that the baryon phase is physical 
at low temperature. As temperature increase, the free energy lines change drastically, see 
Figure |l4|(b). At temperature £t = 0.16, the free energy of baryon phase is lowest in low 
density region. As density increase, there is phase transition between baryon phase and 
quark phase with broken chiral symmetry (first vertical line in |l~4|(b)). After the baryon 



16 



6=0.16 





(c) 



(d) 



Figure 14: Density dependence of free energy for massless quark case with q = 15. To visualize the 
difference of free energy, we plotted F which is defined as F — F — a(T)Q for some a(T). Solid line 
is for baryon phase, red dotted line is for chiral symmetry restored quark phase and dashed purple 
line is for quark phase with chiral symmetry broken. Vertical dotted line denotes phase transition 
point. 



to quark phase transition, there is another phase transition between quark phases: from 
broken to restored chiral symmetry quark phase. At higher temperature, the free energy 
of quark phase is always smaller than that of baryon phase. But there is a phase transition 
from chiral symmetry broken phase to chiral symmetry restored phase at certain density. 



See Figure 14(c), (d). The chiral phase transition between two quark phases coincide with 



the result of previous section, Figure |lfj. Of course, at high enough temperature, the chiral 
symmetry restored quark phase is the only physical phase. These are summarized in the 



phase diagram drawn in Figure [15 . 

This phase diagram has rich structure. At low temperature, baryon phase is always 
physical one in all density region. And chiral symmetry is always broken in this phase. 
However, as we increase temperature, phase transition appears differently depends on den- 
sity. At low density, as temperature increase, baryon phase changes into quark phase but 
chiral symmetry is still broken. If we increase temperature more, then there is chiral sym- 
metry restoration transition in quark phase. On the other hand, at high density, as we 
increase temperature, baryon/quark phase transition and chiral phase transition appears 
at the same time. 

We find that the baryon/quark phase transition temperature decreases as density in- 
crease. However, the decreasing rate is too slow so that phase transition temperature looks 
like constant in the figure. The value of q also affect phase transition. As q increase, both 
phase boundaries move up to larger temperature and larger density region maintaining the 



overall shape. See Figure 15(b). Notice that our phase diagram is very similar to the ones 
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Figure 15: Density dependence of phase transition temperature between quark phase and baryon 
phase (a) with q — 15, (b)with q = 50. 



for Sakai-Sugimoto model studied in ]27|j ; where the baryon phase should be replaced by 
the soliton geometry. Notice that in our model, there is no Hawking Page transition since 
there is no scale other than the temperature and the dilaton contribution is precisely can- 
celed by that of the axion. Our baryon phase is still in the black hole geometry therefore 
the dynamic origin is very different. 

In the case of m q ^ 0, chiral symmetry is explicitly broken so that we might expect 
that there is only baryon/quark phase transition. However, as we discussed in previous 
section, there is another phase transition in quark phase. Look at the short line at the 
up-left region of the phase diagram in Figure 16. It is for m q = 1. 



(t 

1.4 r 



1.2 - 



1.0 - 

o.8 - Quark Phase 

0.6 '- 
0.4 - 
0.2 - 



Baryon Phase 
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Figure 16: Density dependence of phase transition temperature between quark phase and baryon 
phase for m q — 1 and q = 20. 



As we decreases the quark mass, the end point of the phase boundary extends to higher 
density and lower temperature region so that when m q — > 0, we recover the phase diagram 
drawn in Figure [L^. Actually the line becomes long very fast when m q gets to the near zero 
value. Contrarily for q = 0, the line for the transition from black hole embedding (BHE) 
to another BHE exist only in very small density region and it disappear when m q goes to 
zero. Therefore this is very characteristic feature caused by the gluon condensation q. 
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5.2 Relation between chemical potential and density 

Equations of states which are the relations between thermodynamically conjugate variables 
like chemical potential/density quark mass/chiral condensation e.t.c. play important role 
to understand phase structure. In this section, we will discuss the equation of state de- 
scribing the relation between chemical potential and density. Thermodynamically chemical 
potential is defined as derivative of free energy with respect to density 



In quark phase, we can get 

Pquark 



dT, 



quark 



dQ 

i r* &u D7 



27ra' J . dQ 

/CO 
dpdpA t = A t (oo) - A t (p min ). (5.4) 

i-min 

In this case, A t (p m i n ) = because p m in is black hole horizon. It is consistant with usual 
definition of chemical potential(tale of At field) from AdS/CFT correspondence. On the 
other hand, the free energy of baryon phase contains mass of baryon vertex, chemical 
potential should have mass of source. 

"J" baryon 

Hbaryon — qq 

1 [°° 8H m , 1 2 ( ■ 

= A t (oo) - At(p m in) + JT d0n D5- (5.5) 

The last term in ( |5.5|) the baryon mass divided by N c , which is the constituent quark mass. 
Therefore, the chemical potential contains mass of the source. In quark phase, the source 
is the fundamental strings. For the blackhole embedding, D7 brane touches the black hole 
(BH) horizon and the fundamental strings are replaced by deformation of D7 brane. There- 
fore the energy of fundamental strings is contained in that of D7 brane. That's why there 
is no source term, the analogue of the last term in the eq. (|5.5|) ), in the chemical potential of 
the quark phase. The density dependences of chemical potential for m q = embedding are 



drawn in Figure 17. We can see that chemical potentials monotonically increase as density 
increases for any phases. Therefore, there is no thermodynamical instability in m q = 
case. 

The free energy and chemical potential with finite quark mass as functions of density 
are drawn in Figure [j^. Figure fji|(a) shows density dependence of free energy at finite 
temperature with q = 50 and m q = 1. At certain density there is phase transition in quark 
phase. In usual AdS Schwartzschild background without q, there is a thermodynamical 
instability associated with the negative slope branch of the p-Q diagram. But in Figure 
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Figure 17: Density dependence of chemical potential for m q = 0. Red dotted line denotes chemical 
potential for quark phase with chiral symmetry and purple dashed line is for quark phase with broken 
chiral symmethry. The solid line denotes baryon phase. 





(a) 



(b) 



Figure 18: (a) Free energy as a function of density in certain temperature for m q = 1. (b) Density 
dependence of chemical potential, jl = 27ra'/z. In both figure £r = 1.3, q = 50. 



|lq(b), we can see that chemical potential increases at low density and phase transition 
happens before chemical potential begin to decrease. Therefore, chemical potential for 
physical state always increase monotonically as density increase(<9/x/9(3 > 0), hence there 
is no thermodynamical instability. 

As temperature decreases, both quark phase and baryon phase exist and the shape of 
chemical potential is monotonic as function of the density in both phases which is similar 



to Figure 17 
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5.3 Grand canonical ensemble 

In this section, we will discuss the system with grand canonical ensemble. Grand potential 
is thermodynamically defined as 

n = T-nQ. (5.6) 

From definition of 'Hamiltonian density and number density, grand potential for quark 
phase can be written as 

^t-quark — •'quark l^quark^i 



I 



. '"D7\quark phase I " W 

quark phase ' 



(Jp quark phase J J 

I *"D7\ quark phase' ("■') 

The last term in (|5.7| ) is nothing but value of on-shell DBI action for black hole phase. 
It it obvious because 'Hamiltonian density' is defined by Legendre transformation of DBI 
action. Grand potential for baryon phase becomes 

^ ^baryon — •'baryon i^baryorv^i 

= / T~LDl\baryon phase + TT / %D5 _ ( / FQ + TT / ^D5 j 

— / *- D7\baryon phase' \^ - "J 

We can see that energy of source in free energy canceled by one in chemical potential, finial 
form of grand potential is value of DBI action of D7 brane for baryon phase. D5 part does 
not contribute to grand potential. 

In the massless quark case(m,j = 0), the grand potentials of quark and baryon phase 
are drawn in Figure ll^. In each case of the that figure, chemical potential start from non- 



zero value as one can see in Figure 17. It implies that if the chemical potential is not large 
enough to the energy of particle, no particle can be created. Therefore, the system remains 
as vacuum until chemical potential reach to the energy of particle in the system. Therefore 
the value of the chemical potential should be identified as the constituent quark mass. 
Naturally the difference in the chemical potential in baryon phase and that in the realized 
quark phase (between the two quark phases) can be identified as the binding energy. The 
density dependence of the quark mass is plot in figure 19. The phase transition point can 
be identified precisely as the point the binding energy is 0. This point is of course the point 
where baryons melt. The melting point in temperature-density plane is nothing but our 
phase diagram. 



From the figure Figure 19, one can easily read off the critical chemical potential value 
points as one increases the temperature. 

In the case of finite quark mass, chiral symmetry is always broken similar to the canon- 
ical ensemble. We have only two phases, the baryon phase and quark phase. However, we 
find that the density dependence of chemical potential and the chemical potential depen- 
dence of grand potential behave similarly to the massless case except the absence of chiral 
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Figure 19: Chemical potential dependence of grand potential(m 9 = and q = 15). Red dotted 
line denotes to grand potential for chiral symmetry restored quark phase purple dashed line to 
grand potential for quark phase without chiral symmetry and solid line to baryon phase 
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Figure 20: Phase boundary in grand canonical ensemble in the case of (a) m q = 0, (b) m q = 1. 



phase transition. The phase structure based on above discussion is drawn in Figure 20. 
Similar to phase diagram in canonical ensemble, the temperature of phase transition be- 
tween quark and baryon phase is almost constant. We also find the transition in quark 
phase in small density region, but this line is not shown in this diagram. 



5.4 Chiral condensation 

One of the most important observable in QCD is the chiral condensation (CC) and its 
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density dependence. It measures the dynamical mechanism for creating the mass of the 



hadrons. Figure 21 shows the density dependence of the chiral condensation. In the baryon 
phase, chiral symmetry is always broken and there is a quark phase where chiral symmetry 
is broken. Our result shows that CC is increasing in baryon phase and decreasing in the 
relevant quark phase. 
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Figure 21: Density dependence of the value of chiral condensation < qq > for massless quark case 
with q — 15. Thick line indicates physical phase from free energy and there is phase transition 
along vertical dashed line. 



6. Summary and Discussion 

In this paper, we study phase structure of a holographic QCD based on D3/D-instanton 
background. The phase transition here is for confinement/deconfiiement of quark not the 
gluons. While the Hawking Page transition in usual AdS/CFT correspondence describes 
the dynamics of the gluon, confinement/deconfinement phase transition of quarks are de- 
termined from the interaction between the geometry and the compact brane dynamics. 
Namely it is question of existence of the baryon vertex. The D3/D-instanton background 
have quasi-confining nature as defined in the introduction. In zero quark mass case, there 
is a chiral transition in quark phase. The number of D-instanton, identified as the expec- 
tation value of gluon condensation played, essential role in breaking the chiral symmetry. 
The phase boundary is depends on value of q. 

A few remarks are in order: 

The first one is on the assumption that there is no Hawking-Page transition. Main 
reason is that the geometric transition should be discussed in the Einstein frame, where 
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dilaton action is cancelled by that of axion therefore we can not find any effect of the 
dilaton condensation scale in bulk free energy level. Therefore it is completely the same as 
the D3 brane without D-instanton. However all the probe D-brane dynamics is affected by 
the presence of the dilaton factor. That is why we get non-trivial result. This is interesting 
and subtle point and in the main sections we just assumed that there is no Hawking-Page 
transition. 

Second one is about the role of the Chern-Simons term. In ref. |J, the authors intro- 
duced Chern-Simons term such that it cancel out the dilaton effect of the brane embedding 
dynamics. However, we could not find solution with such behavior. We found that the 
Chern-Simons term is a total derivative so that it can contribute to the charge but not to 
the equation of motion. This allows us the non-trivial effect of the gluon condensation on 
the embedding of D7 brane as well as on that of D5. 

The third one is that due to the probe nature of the embedding dynamics, our cal- 
culation is not trusted in extreme high density regime. One need to take care of the 
backreaction of the geometry. 

The 4th is the Euclidean nature and duality in such background. For the Euclidean 
configuration, there is no state/geometry correspondence. However, gauge gravity duality 
is still there. The correspondence of D-instanton in type IIB and and Yang-Mill theory 



instanton was discussed in [28] as well as many other papers. The relation of AdS/CFT and 
multi-instanton gauge theory is considered in well known works of Dorey et.al. 129(1 . So the 
gauge theory dual of the uniform distribution of such D-instanton is not hard to imagine 
and it should be the Yang-Mill theory with uniform distribution of instanton charge. 

The 5th is the Infra-red singularity. The background has the IR singularity since the 
dilaton factor has log singularity at the horizon. Usually one need to prevent any D-brane 
probe to approach such singular region. Also this makes the bulk action divergent and 
we would need IR cut off. However, in our background, both the bulk gravity action as 
well as the DBI action are regular at the horizon. There is no divergence. There are 
reasons for this, i) The finite temperature version has a regular horizon and the essential 
singularity is hidden in the horizon, ii) The finiteness of the bulk action is partly due to 
the supersymmetric construction of the action and the ansatz eq.(3.1) by which the action 
of dilaton is cancelled precisely by that of axion. iii)The finiteness of the probe brane 
action is due to detailed structure of DBI action. Looking at eq. (4.7), the log divergence 
of dilaton factor e* is cancelled by the presence of the w_ factor, which has a zero. The 
thermodynamics is based on the calculation of the DBI action for the actual configuration 
of the Speaking more physically, the background act net repulsive force to the D-brane 
which perform the dynamic censorship and gives finite values of DBI action. 

In the background with q = 0, the meson spectral function does not have interesting 
feature. But due to the repulsive nature of the force acting on the probe brane, there is 
some reason to expect non-trivial spectral behavior in this case. It would be very interesting 
to calculate the meson spectrum in the chiral symmetry broken quark phase. 
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A. Chern-Simons term 

We will put probe D7 brane, then 8-form field from hodge dual of axion can interact 
with D7 brane world volume. We check it in zero temperature limit and extend to finite 
temperature case. 



In the limit T — y 0, (2^) becomes 

ds = e*/ 2 { ^ (-dt 2 + d&) + ^ (dr 2 + r*ctt§) } 

e* = l + ^, X = -e-* + Xoo- (A.l) 

We are interested in interaction between dual field of axion and probe D7 brane, we change 
the metric into the direction along and perpendicular to D7 brane, 

ds = e */ 2 | ^ (-dt 2 + dx 2 ) + ^ {dp 2 + p 2 dQ 2 + dy 2 + y 2 # 2 ) j , (A.2) 

where r 2 = p 2 + y 2 . To get dual field, we introduce vielbein 

e * = e */4ldt, e* = e*/ 4 -^, e P = e*/ A -dp, 
R R r 

e ^a = e^Xpdn^ e y = e */^ d ^ e * = e */4* d0 (A3) 

/y> fy* ty* 

where tilde denotes to flat index. The field strength for axion field is 

= e ~7T d P + e ~7T d y 
op ay 

dpR + 8yR ( } 

Hodge dual of this field strength is 

(9} 9! dp R txiCi 3 y4> 

1 - 5 $/4^ r j j ^ p o 3 A J 
9! <9y R txia 3 p<i> 
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e — — p ydt A dx A df£ 3 A dy A d<; 
dp 

-e —-p ydt A dx A d£l 3 A dp A < 
oy 

— — p 3 ydt A dx A d^ 3 A dy A d</> 
dp 

dy 



p 3 ydt A df A df) 3 A dp A • 



where we are using the convention 



"txipQ-iycp 



+ 1. 



By substituting e*, we get 



\9) 



(p 2 + y 2 ) 3 

4gp 3 y 2 



(p 2 + y 2 ) 3 

Now, we want to get 8-form potential such that 



dt A dx A dil 3 A dy A d<f> 
dt A dx A d0 3 A dp A dip. 



*(9) 



dC, 



(»)■ 



Assuming 



C, 



(8) 



/(p, y, 0)dt A dx A d$7 3 A dy 
+y(p, y, 4>)dt A dx A dQ, 3 A dp, 



we get 



df 
dC (8 ) = -^-dt A dx A dn 3 Ady Ad(p 
v ' o<p 

df 
—Trdt A dx A dtt 3 A dp Ady 
op 

+Trrdt A dx A dQ, 3 A dp A d(j) 

oq> 

da 
+ Trdt A dx A d0 3 A dp A dy. 

oy 



(A.5) 



(A.6) 



dg 


(p 2 + y 2 ) 3 

4qp 3 y 2 


d(j> 
dp 


(p 2 + y 2 ) 3 ' 
oy 



By integrating / and g with respect to <f), we get 

4gp 4 y 



f(p,y,4>) 



(p 2 + y 2 ) 



2\3 



+ K 



(A.7) 



(Ai 



(A.9) 



(A.10) 



By comparing ( A.7 ), we get the condition for f(p,y,(j>) and g(p,y,(f>) as follows; 

df ^qpS 



(A.ll) 
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g(p,y,<i>) 



±qp y 



3„2 



but from the last condition of ( A.ll 
be written 



( p 2 + y 2)3Vr 

we get 0i = <f> 2 



+ 2 ), (A.12) 

= (fro. Finally, 8- form potential can 



a 



w 



_ 4qp y 

{ P 2 + y 2 Y 
4gp 3 y 2 , 

(^2 + 1,2)3 W 



'0 + <fro)dt AdfA dil 3 A dy 
+ 0o)dt AdxA df2 3 A dp. 



(A.13) 



If we fix location of D7 brane along cfr direction at 
derivative i.e. 



a 



(8) 



4gp 4 y ( 

(p 2 + y 2 ) 3 
4gp 3 y 2 , , 

(p2 + y 2)3^ 

(0* + O ) d 



= 0*, then ( |A.13 ) becomes total 
* + (fro)dt AdiA dfi 3 A dy 



+ 4>o)dt AdfA d(7 3 A dp 

„4 

rdt Adz AdS7 3 



<7P 



(A.14) 



V + y 2 ) 2 

so that C(g) is a pure gauge whose field strength is zero. Furthermore we can always choose 
00 = ~ < fr* + 2vr, then the Chen-Simons term in ( |3.3j ) becomes 

„4 

-dt A dx AdS7 3 



S CS = P7 j d 8 a(2TT)d 



qp 



(27r)fi 7 v 4 n 3 



qp 



(p 2 + y 2 ) 2 
,i 



(/o 2 + y 2 ) 



2^2 



p=oo 



g(27r) j U7T40 3 



2^(-D- 



(A.15) 



It is nothing but the half of D-instanton number calculated in §. This is completely 
satisfactory since D7 brane world volume can captures the flux of the upper hemi-sphere. 
Since the Chern-Simons action is locally a total derivative term, it does not contribute to 
the equation of motion. 

B. Hamiltonian density of baryon vertex 

We start from the action for D5 brane with Chern-Simons term (|4.7l), 

Sd5 = Sdbi + Scs 

= -p 5 I e~* v / -det( 5 + 2vtq'F) + p 5 f A {1) /\G 



(5) 



r 5 / dtd9 sin 4 6 



,* 



UJ- 
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dtd6£ D5 , (B.l) 

where 

T 5 = fi 5 n A R 4 r T , F = 2ira'F te , A t = 27ra'A t . (B.2) 

We denote £^5 to Lagrangian density. The displacement can be obtained by derivative of 
Lagrangian density with respect to A' t , 

sin 4 0F . x 

= -2toV 5 — . (B.3) 



Then the equation of motion for gauge filed can be written as 

d e D(6) = -4 sin 4 6, (B.4) 

where D(9) = 2n a'r ■ This equation plays a constraint in the action, then the action of D5 
brane becomes, 

S D 5 = S DB i + t 5 dtd04 sin* 6 A t 



Sdbi - t 5 I dtd6{d e D{9))A t 

Sdbi-t 5 I dtdODF, (B.5) 



where we take integration by part in last procedure. This is nothing but the Legendre 
transformation of DBI action of D5 brane. After substituting ( |B.3| ) in the the action, we 



define 'Hamiltonian density' (4. 9) as follows 



S D5 = -T 5 J "*Wy^ (e 2 + e)\/D(0) 2 + sm 8 e 

dtd6U D5 . (B.6) 



The definition of Hamiltonian density is consistent with one of probe D7 brane and inte- 
gration of this with on-shell solution gives free energy of D5 brane. Next, by integrating 
( |B.4| ) , we get 

D(0) = !(i^_#) + jj s in0cos6< + sin 3 0cos#, (B.7) 

where the integration constant v determines the number of fundamental strings on each 
pole, i. e. vN c strings are attached to the south pole and (1 — v)N c strings to north pole 
of D5 brane. Here, we assume that all the fundamental strings are attached on the north 
pole, we set v = 0. 
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C. Force balance condition 



In this section, we derive force balance condition ( 4,11 ). In our mode, the end points 
of fundamental strings play role of source of U(l) gauge field. Due to the tension of 
fundamental string, there can exist cusps on probe brane world volume. By calculating 
force at the cusp of each brane, we can estimate the behavior of probe branes. 

The force at the cusp can be obtained by taking small variation of 'on-shelP free energy 
with respect to U c 

"./"on—shell 



SU r . 



(C.l) 



where U c is position of cusp and the free energy is integration of Hamiltonian density of 
probe brane 

T = I dpHonshelb (C.2) 



The Hamiltonian density is a function of U and U' , we can write the variation as follows; 

^on-shcll = o"H(U,U ;/?) on -shell 



Cr£on— shell r- rT , Cr£on— shell rTT i 

-oil H tttt: oil 



dU 

5%on-shell 



d 



d_ 

dp 
d_ 
dp 



au 6u+ d P 

^Hpn-shcll 
8U' 

^^on-shell 
8U' 



au' 

dU 



on— shell 



au' 



5U 



d ( ^Hpn-shell 

dp V dU> 



Finally, the force at the cusp is 



F r 



5U 
5U 

dp 
dp 



+ 



d'Hon-shell d_ ( ^on-shell 

dU ~ dp\ dU' 



5U 
8U 



(C.3) 



sn 



on— shell 



SU C 



dp 



dTirm- 



on— shell 



su 



au> su r 



<9%on-shell 



au> 



(C.4) 



U=Uc 



The force at the cusp of single D5 brane can be calculated from (C.4), 



F m = N c T F1 i 



'e*wl 



2 uj. 



Ve+l 



a 



(C.5) 



e=Sc 



where Tp\ is tension of fundamental string \/2ira' and ^ c is position of cusp of D5 brane. 
With a same manner, we can get force at the cusp of probe D7 brane, 



Fd7 = Tpi\ 



l e*u?_ 



Qy 



2 W+ y/l + y2 



y = yc, 



(C.6) 
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where y c is position of D7 brane at the cusp(p = 0) and Q is total number of U(l) source. 
The force at the cusp of D5 and D7 is always smaller than the fore of fundamental strings. 
Then strings pull each brane until the length of string becomes zero and force between 
two branes is balanced. We can get condition of slope of D7 brane from the force balance 
condition. 

= F D7 {Q) + N B F D5 
= F D7 {Q) + ^-F m 



f. (C.7) 



30 



References 

[1] E. Witten, Adv. Theor. Math. Phys. 2, 505 (1998) [arXiv:hep-th/9803131]. 

[2] E. Witten, JHEP 9807, 006 (1998) [arXiv:hcp-th/9805112]. 

[3] C. G. . Callan, A. Guijosa, K. G. Savvidy and O. Tafjord, Nucl. Phys. B 555, 183 (1999) 
[arXiv:hep-th/9902197]. 

[4] Y. Seo and S. J. Sin, JHEP 0804, 010 (2008) [arXiv:0802.0568 [hcp-th]]. 

[5] T. Schafcr and E. V. Shuryak, Rev. Mod. Phys. 70, 323 (1998) [hcp-ph/9610451]. 

[6] H. Liu and A. A. Tseytlin, Nucl. Phys. B 553, 231 (1999) [arXiv:hep-th/9903091] . 

[7] S. S. Gubser, arXiv:hcp-th/9902155. 

[8] Y. Kim, B. H. Lee, C. Park and S. J. Sin, Phys. Rev. D 80, 105016 (2009) [arXiv:0808.1143 
[hep-th]]. 

[9] K. Ghoroku, T. Sakaguchi, N. Uckusa and M. Yahiro, Phys. Rev. D 71, 106002 (2005) 
[arXiv:hep-th/0502088]. 

[10] K. Ghoroku and M. Ishihara, Phys. Rev. D 77, 086003 (2008) [arXiv:0801.4216 [hcp-th]]. 

[11] S. J. Sin and Y. Zhou, JHEP 0905, 044 (2009) [arXiv:0904.4249 [hcp-th]]. 

[12] K. Y. Kim, S. J. Sin and I. Zahed, "Dense hadronic matter in holographic QCD," 
[arXiv:hcp-th/0608046] . 

[13] N. Horigomc and Y. Tanii, JHEP 0701, 072 (2007) [arXiv:hep-th/0608198]. 

[14] S. Nakamura, Y. Seo, S. J. Sin and K. P. Yogcndran, J. Korean Phys. Soc. 52, 1734 (2008) 
[arXiv:hep-th/0611021]; 

S. Nakamura, Y. Seo, S. J. Sin and K. P. Yogcndran, Prog. Theor. Phys. 120, 51 (2008) 
[arXiv:0708.2818 [hcp-th]]; 

[15] S. Kobayashi, D. Matcos, S. Matsuura, R. C. Myers and R. M. Thomson, JHEP 0702, 016 
(2007) [arXiv:hep-th/0611099]. 

D. Matcos, S. Matsuura, R. C. Myers and R. M. Thomson, JHEP 0711, 085 (2007) 
[arXiv:0709.1225 [hcp-th]]. 

[16] G. W. Gibbons, M. B. Green and M. J. Perry, Phys. Lett. B 370, 37 (1996) 
[arXiv:hep-th/9511080]. 

[17] A. Kehagias and K. Sfctsos, Phys. Lett. B 456, 22 (1999) [arXiv:hcp-th/9903109]. 

[18] K. Ghoroku and M. Yahiro, Phys. Lett. B 604, 235 (2004) [arXiv:hcp-th/0408040]. 

[19] M. Kruczcnski, D. Mateos, R. C. Myers and D. J. Winters, JHEP 0307, 049 (2003) 
[arXiv:hep-th/0304032] . 

[20] J. Babington, J. Erdmenger, N. J. Evans, Z. Guralnik, I. Kirsch, Phys. Rev. D69, 066007 
(2004). [hcp-th/0306018]. 

[21] M. Kruczenski, D. Mateos, R. C. Myers and D. J. Winters, JHEP 0405, 041 (2004) 
[arXiv:hcp-th/0311270]. 

[22] J. Erdmenger, A. Gorsky, P. N. Kopnin, A. Krikun and A. V. Zayakin, JHEP 1103, 044 
(2011) [arXiv: 1101. 1586 [hep-th]]. 



31 



[23] M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B 147, 385 (1979). 

[24] D. Matcos, R. C. Myers and R. M. Thomson, Phys. Rev. Lett. 97, 091601 (2006) 
[arXiv:hep-th/0605046] . 

D. Matcos, R. C. Myers and R. M. Thomson, JHEP 0705 (2007) 067 
[arXiv:hep-th/0701132]. 

[25] N. Evans, A. Gebauer, K. Y. Kim and M. Magou, JHEP 1003 (2010) 132 [arXiv:1002.1885 
[hep-th]]. 

[26] O. Bergman, G. Lifschytz and M. Lippert, "Holographic Nuclear Physics," 
[arXiv:0708.0326]. 

[27] O. Aharony, J. Sonnenschein and S. Yankielowicz, Annals Phys. 322, 1420 (2007) 
[arXiv:hep-th/0604161]. 

[28] M. Bianchi, M. B. Green, S. Kovacs and G. Rossi, JHEP 9808, 013 (1998) [hep-th/9807033]. 

[29] N. Dorey, T. J. Hollowood, V. V. Khozc, M. P. Mattis and S. Vandoren, Nucl. Phys. B 552, 
88 (1999) [hep-th/9901128]. 



32 



